Recoil effects of photoelectrons in a solid 
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ABSTRACT 

High energy resolution C Is photoelectron spectra of graphite were measured at the excitation 
energy of 340, 870, 5950 and 7940eV using synchrotron radiation. On increasing the excitation 
energy, i.e., increasing kinetic energy of the photoelectron, the bulk origin C Is peak position shifts 
to higher binding energies. This systematic shift is due to the kinetic energy loss of the high-energy 
photoelectron by kicking the atom, and is clear evidence of the recoil effect in photoelectron emission. 
It is also observed that the asymmetric broadening increases for the higher energy photoelectrons. 
All these recoil effects can be quantified in the same manner as the Mossbauer effect for 7-ray 
emission from nuclei embedded in crystals. 

PACS numbers: 79.60.-i, 79.20.-m, 79.20.-m 



Photoelectron spectroscopy is widely used for the 
study of electronic structure of solids [lj. The binding 
energy Eb of the electron is often calculated from the 
photoelectron kinetic energy using the following equa- 
tion: 
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where Ekm is the measured electron kinetic energy, hv 
is the photon energy for excitation, and <j) is the work 
function. This procedure overlooked recoil effects, which 
lead to part of the kinetic energy being imparted to the 
emitting atom. As a result, the binding energy deter- 
mined in this way is greater than the true binding en- 
ergy. This effect is very small for vacuum ultraviolet and 
soft x-ray photoelectron spectra, so that recoil effects can 
be safely neglected. For photoelectrons with 1000 eV ki- 
netic energy emitted from carbon atom, the recoil energy 
is estimated to be only ~ 45 meV, although Kukk et al. 
succeeded recently to observe small deviation in spectral 
shape of vibronic lines in gaseous methane, which have 
been attributed to a recoil effect Q. 

The momentum transfer at recoil is a fundamental pro- 
cess observed in experiments of neutron and x-ray scat- 
tering high-energy electron backscattering etc. 
For photoelectron emission, Domcke and Cederbaum Q 
predicted that the recoil effect can be observed as a spec- 
tral modification for gaseous molecules with light atoms. 
Quite recently, Fujikawa et al. Q evaluated the amount 
of shift and broadening of core-level photoelectron spec- 
tra, as well as for electron backscattering, due to recoil 
effect in solids. It is noted that at keV energies, since the 
momentum of an electron is much larger than that of a 
photon of the same energy, and the transferred momen- 
tum is largely that of the emitted electron, it should be 



possible to detect recoil effects in photoelectron emission 
with keV energies. 

In the last few years, hard x-ray photoelectron spec- 
troscopy with the excitation energy of 6-8 keV has 
been realized using high brilliance synchrotron radia- 
tion [H, 0, [n| , resulting in useful studies on semiconduc- 
tors and correlated materials [ll|, [12] • Since the achieved 
energy resolution is quite good (AE < 80 meV), it gives 
us an opportunity to investigate recoil effects in a solid. 

In this study, we measured the C Is core level pho- 
toelectron spectra of highly oriented pyrolytic graphite 
(HOPG), using high energy resolution (AE = 100 ~ 120 
meV), with the aim of investigating recoil effects due to 
photoelectron emission in a solid. HOPG is a suitable 
material for study of the recoil effect because of the light 
atomic mass of carbon. Furthermore, reliable studies of 
surface-bulk core level shift of HOPG have already been 
established recently [l3|, 14, 15, providing a very 
suitable reference for this study. The experimentally ob- 
tained spectra in the present study confirm the surface 
core level shift. For the bulk derived feature, the spectra 
exhibit systematic kinetic energy loss and the anomalous 
asymmetric broadening on increasing the excitation en- 
ergy. All these features are clear evidence of recoil effect, 
and can be theoretically quantified in the same picture 
as the Mossbauer effect for 7-ray emission. 

Measurements of C Is photoelectron spectra of HOPG 
were performed at SPring-8 using synchrotron radiation. 
Hard x-ray photoelectron spectra at the excitation en- 
ergy of 5950 and 7940 eV were measured at the undulator 
beamline BL29XU using a hemispherical electron energy 
analyzer, SCIENTA R4000-10kV. Details of the appara- 
tus including x-ray optics are described in refs. [9l.llfl.ll8l]. 
Soft x-ray spectra at the excitation energy of 340, 870 eV 
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were measured at the undulator beamline BL17SU [l9[ 
using an electron analyzer, SCIENTA SES-2002. Clean 
surfaces of HOPG were prepared by peeling off an ad- 
hesive tape at a pressure of 10~ 8 Pa for all measure- 
ments. All the measurements on graphite reported here 
were carried out at room temperature. The energy scale 
of the spectra were calibrated very accurately (<10 meV) 
by measurements of the Au 4f core levels and the Fermi 
edge at room temperature. The total instrumental en- 
ergy resolution (AE) for the soft x-ray and hard x-ray 
spectra were determined as 100 meV and 120 meV by 
fitting the Fermi-edge profiles of Au measured at 20 K. 

C Is core level spectra of HOPG measured at room 
temperature with soft x-ray excitation (340 and 870 eV) 
and hard x-ray excitation (5950 and 7940eV) are shown 
in Fig. 1(a). The spectra are obtained effectively in 
normal emission geometry. The photoelectron detection 
angle relative to the sample surface for the soft x-ray and 
hard x-ray spectra were 90 ° and 85 ° , respectively. The 
peak position of the spectrum obtained with hv=870 eV 
shifts to lower binding energy (284.39 eV)in comparison 
with that of his=M0 eV spectrum (284.47 eV). This shift 
is due to the difference in the probing depth of photoelec- 
trons between these excitation energies. With increase of 
the kinetic energy, the probing depth of a photoelectron 
becomes larger. The spectrum of hi/=340 and hv=870 
eV are dominated by the surface and bulk components, 
respectively. The observed surface-bulk core level shift is 
consistent with the results reported by Balasubramanian 
et al. [IH , using photon energies of 300-348 eV, and a to- 
tal energy resolution of about 50 meV. It was shown that 
surface-bulk core level splitting results in a weak bulk fea- 
ture at a lower binding energy compared to the surface 
derived feature. In the present study, our /i^=340 eV 
spectrum is similar to that reported by Balasubramanian 
et al. [TBI, while the ft,z^=870 eV spectrum is dominated 
by the bulk derived feature. 

With increase of the excitation energy, i.e., the ki- 
netic energy of the C Is photoelectron, the peak position 
shifts to higher binding energy side. It is also obvious 
that asymmetric broadening becomes much wider for the 
higher energy photoelectrons, while the total instrumen- 
tal energy resolution for these spectra is almost the same. 
For the possible origin of these spectral changes, the ef- 
fect of associated elementary excitation such as plasmons 
is discarded, since the peak position itself shifts, depend- 
ing on the kinetic energy. The asymmetric line shape of 
the C Is core level spectra obtained with soft x-ray ex- 
citation has been discussed in relation with semimetallic 
character of gra phit e, and can be fitted by the Doniach- 
Sunjic function [l3l . Q, [H, This possibility can be 
also excluded for the same reason. Thus, we are led to a 
picture based on recoil effects for explaining the observed 
peak shift and broadening, which depends on the kinetic 
energy of the photoelectron. 

For an atom with mass M in free space, the recoil 
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FIG. 1: (color online) (a) Photon energy dependence of C Is 
core level spectra of graphite. The soft x-ray (hv = 340, 870 
eV) and hard x-ray (hv = 5950, 7940 eV) are measured at the 
emission angle of 90 ° and 85 ° relative to the sample surface, 
(b) Theoretically obtained spectra taking into account the 
recoil effect in a Debye model with huji, t n —75 meV. 

energy SE is simply estimated from the momentum con- 
servation as SE = (m/M)Ekin, where m is the electron 
mass. For a carbon atom with the mass ratio m/M = 
1/22000, SE becomes as large as 0.36cV for E km =8keV. 
In the solid, this recoil energy is absorbed by the phonon 
bath, resulting in the excitation of phonons. The zero- 
phonon transition corresponds to the event in which the 
recoil energy is transferred to the center of mass motion 
of the total crystal. This is essentially the same as the 
Mossbauer effect in the 7-ray emission from nuclei em- 
bedded in crystals [iol ]. 

Without loss of generality, we assume that a core elec- 
tron of the carbon atom located at the lattice point R° 
is ejected by the x-ray irradiation. The actual nuclear 
position R may deviate from its equilibrium value R° as 
R = R° + u because of the thermal and zero point vi- 
brations. The interaction Hamiltonian with the x-ray is 
given by Hj = (a + a^e-p, aside from irrelevant factors, 
where a is the annihilation operator for the x-ray with en- 
ergy hv and the polarization vector e. The momentum of 
photon is neglected here, since it is an order of magnitude 
smaller than the that of the emitted electron in this x- 
ray energy region. The initial state of transition is given 
by \^>i) = \hv) (g> \4> c ) ® where \hv) is the one-photon 
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state, \i(j c ) is the core electron state with energy e c , and 
\i) is a phonon state of the crystal. The crucial point of 
the theory is that the wave function of the core electron 
is given by the form (r \ip c ) = 4 } c(r — R). This is the adia- 
batic approximation, and the recoil effect results directly 
from this functional form. On the other hand, the final 
state of the transition is given by I*/) = |0) ® \tpk) ® |/), 
in which |0) is the vacuum of the photon, \ipk) is the plane 
wave of the electron, (r \ipk) = (2-/t)~ 3 / 2 cxp(ik ■ f), with 
energy H 2 k 2 /2m , and |/) is a phonon state. 

Using the functional forms for and |\&/), and 

changing the integration variable from r to r — R , we 
find 



<*/|ff/|*i> 



-ik-R 



in which 



ft = (2tt)- 3 / 2 J d 3 re- lk p ( j ^ c (r). 

The above expression of the interaction Hamiltonian co- 
incides with that for the Mossbauer effect [2l|, so that 
we readily obtain the expression for the photoelectron 
spectrum as a function of the relative binding energy E 
measured from the recoilless value, 
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in which the generating function F(t) is given by the 
canonical average, 



Fit) = ( e ik -' a W e ~ ik - s 



with u(t) being the Hciscnberg representation of u at time 
t, and r is the lifetime broadening factor. For harmonic 
crystals, F(t) can be written in the closed form, F(t) = 
exp[G(t)] where 

G(t) = JI a 9 {( 2n ( w <z) + 1 ) ( C0SUJ q t - 1) - isincjqt} , 



with 
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and 
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in which g is the abbreviation for the wave vector and 
the branch index of phonons, M is the mass of unit cell, 
N is the number of unit cells, and ff q is the polariza- 
tion vector of the phonon. The photoelectron spectrum 
1(E) is the convolution of a structure function with the 
Lorentz function. The former consists of the recoilless 



line (the zero-phonon line) and its phonon side bands. 
Unlike the Mossbauer effect, one can control the so-called 
Debye- Waller factor continuously from almost recoilless 
to strong coupling by changing the energy of the x-ray. 
The spectrum may depend also on the relative angle 
of the k vector of the emitted electron because of the 
anisotropy of the phonon spectrum. 

In actual calculation, we adopt an anisotropic Debye 
model for graphite. The stretching (in-plane) mode and 
the bending (out-of-plane) mode are assumed to be in- 
dependent. The surfaces of constant frequency for each 
mode have an prolate spheroidal form in q space reflecting 
the highly anisotropic dispersion relation for the in-plane 
propagation and the out-of-planc propagation (2^]. Then, 
G(t) is given by 



G(t) 



duj(e lut - 1){ J„(w) cos 2 9 + J h (uj) sin 2 6} 



(1) 

for the emission angle 9 relative to the surface, where 
J\(u>) (A = s, b) are given by 

Jx(f») = j-{(n(u;) + l)D x (cj)e(Lj)+n(\u\)D x (\L;\)e(-Lu)} 

with the step function Q(lj). Here, D\(lo) is the density 
of state per atom, 



Dx(uj) 
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(2) 



Su} \,d u \,c i otherwise. 



The density of state has a characteristic of three- 
dimensional Debye model in low frequency (to < lu\,c), 
but of two-dimensional one at high frequency, where 
only the in-plane propagating modes contributes to 
D\(w). The Debye cut-off frequencies for the stretch- 
ing mode fno S} D and the bending mode huih,D are esti- 
mated from the frequency at if-point of LA mode and 
ZA mode, respectively. From the experimental disper- 
sion curves [23j|. we fix, hw s _D =150mcV, hujf, d =75meV, 
hw s ,c =6.3meV, fuvb.c =12.5mcV. 

The theoretical spectra of the C Is normal emission 
photoelectrons of graphite were calculated and are plot- 
ted as a function of the recoil energy in Fig. 1(b). The 
lifetime broadening (full- width at half maximum) is taken 
to be 160 meV as is known for graphite [3, EH- The 
spectra are then convoluted with a Gaussian function 
corresponding to the experimental resolution of 120 meV. 
Theoretical spectra reproduce the experimental spectra 
fairly well without any adjustable parameters. Note that, 
at 870 eV excitation, the peak shift is quite small, but 
phonon excitations cause the asymmetric broadening. 

Fig. 2 (a) shows the experimental C Is core level spec- 
tra of HOPG measured with 7940 eV excitation at the 
emission angle of 85 ° (the same as that in Fig. 1(a)) 
and 30 ° relative to the sample surface. The peak slightly 
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FIG. 2: (color online) (a) C Is core level spectra of graphite 
measured at the photon energy of 7940 eV at the emission an- 
gle of 85 ° (normal) and 30 ° (grazing) relative to the sample 
surface, (b) Theoretically obtained spectra with the Debye 
energy fkJb,D=75 meV and tkj a ,D = i-50 meV for the bending 
and the stretching mode, respectively. 

shifts to lower binding energy and becomes broader for 
the grazing emission spectrum. Theoretical spectra are 
shown in Fig. 2 (b). The theoretical spectra reproduce 
well the observed emission angle dependence. The larger 
spectral width in the grazing angle emission than the 
normal one is attributed essentially to the large Debye 
cut-off frequency uj s .d f° r the stretching mode, which is 
roughly twice that for the bending mode u)t,D- 

It is remarkable that the gross spectral features in Fig. 
1 can be understood by a semiclassical approximation [7| • 
We have checked that the spectra can be approximately 
fitted by a Gaussian line shape, using an atomic value for 
the shift SE, and a width given by the second moment 
due to the Doppler broadening of 25EksT. However, 
the semiclassical picture fails in reproducing the asym- 
metries of the observed line shapes, which are signatures 
of the quantum nature of phonons. The anisotropic line 
shape shown in Fig. 2 also indicates the necessity of a 
full quantum mechanical analysis taking into account the 
solid state effect for such a high resolution spectroscopy. 
Furthermore, the good agreement between the experi- 
mental data and the theoretical calculation tells us that 



the multiple-scattering effect is small. 

In conclusion, recoil effects in photoelectron emission 
from solid were observed in the high energy C Is core 
level spectra of graphite. The observed spectral shapes 
are successfully quantified in the same manner as the 
Mossbaucr effect for 7-ray emission from nuclei embed- 
ded in crystals. The recoil effect exists, and is measur- 
able, making it an important aspect of high-energy pho- 
toelectron spectroscopy. 
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